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Abstract —In this paper, a quantum filter for estimating the 
states of a non-Markovian quhit system is presented in an aug¬ 
mented Markovian system framework including both the quhit 
system of interest and multi-ancillary systems for representing the 
internal modes of the non-Markovian environment. The colored 
noise generated hy the multi-ancillary systems disturbs the qubit 
system via a direct interaction. The resulting non-Markovian 
dynamics of the qubit is determined by a memory kernel function 
arising from the dynamics of the ancillary system. In principle, 
colored noise with arbitrary power spectrum can be generated by 
a combination of Lorentzian noises. Hence, the quantum filter can 
be constructed for the qubit disturbed by arbitrary colored noise 
and the conditional state of the qubit system can be obtained by 
tracing out the multi-ancillary systems. An illustrative example is 
given to show the non-Markovian dynamics of the qubit system 
with Lorentzian noise. 

I. Introduction 

The qubit is a fundamental unit of quantum computation 
and information, which has been a research focus for the past 
two decades III. The nature of quantum superposition endows 
a qubit with the capability of carrying more information than 
a classical bit. Hence, qubit-based quantum computation can 
speed up calculations in a suitable algorithm. 

Many potential systems for constructing qubits have been 
investigated, e.g., nuclear magnetic resonant systems, super¬ 
conducting systems and quantum dots systems III. In these 
potential systems, solid-state systems have been paid more 
attention due to long coherence times, scalability, and conve¬ 
nient operations and readouts 121 ■ However, due to the memory 
effects in the nature of solid-state systems, non-Markovian 
dynamics of qubit systems have to be dealt with 13, where 
the commutation relation of the colored noise is determined 
by a memory kernel function of the environment a, HI. 
Correspondingly, when considering the non-Markovian effects 
of colored noise in classical control engineering, a whitening 
filter is often used by appending the state of the colored noise 
model to that of the plant, resulting in a Markovian dynamics 
of the augmented states driven by white noise |6l, iTj- 
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In this paper, we represent a non-Markovian qubit system 
in an extended Markovian system framework. In particular, 
n ancillary systems defined on a Hilbert space i}®” are 
introduced to play the role of the internal modes of the 
non-Markovian environment converting white noise to colored 
noise. The structure of the ancillary systems determines the 
spectrum of the colored noise. Supposing the qubit system is 
defined on a Hilbert space t) and the noise field is defined on 
a Fock space the Markovian evolution of the augmented 
system is defined on t) ig) 1}®" (g) Such an approach was 
named as a pseudo-mode model for non-Markovian quantum 
systems HI, ii and was applied to model energy transfer 
process in photosynthetic complexes nni. Similarly, the dy¬ 
namics of non-Markovian quantum systems can be described 
by a hierarchy equation approach HU where parts of the 
equations describe the pseudo-mode dynamics. This has been 
applied to the indirect continuous measurement of a non- 
Markovian quantum system lfT2l . However, this pseudo-mode 
approach has not yet been systematically described so as to be 
compatible with quantum control theory, e.g., quantum filtering 
theory. 

The multi-ancillary systems for representing the internal 
modes of the environment are described by quantum stochastic 
differential equations (QSDE) in this paper. For the fictitious 
output of each ancillary system with a Lorentzian power 
spectrum, the spectrum of the colored noise arising form the 
ancillary systems is a combination of Lorentzian ones such 
that colored noise with an arbitrary power spectrum can be 
approximately generated ns. This colored noise disturbs the 
qubit system via their direct interactions such that the dynamics 
of the qubit system can be described by a quantum stochastic 
integral differential equation (QSIDL) or a non-Markovian 
Langevin equation. In addition, the augmented model of the 
non-Markovian quantum system can be conveniently described 
by an {S, L, H) description in an extended Hilbert space which 
is compatible with quantum filtering theory. By applying a 
probing field into the qubit system, a quantum filter for the 
non-Markovian qubit system can be constructed. Due to the 
output quadrature satisfying a non-demolition condition, the 
augmented system state can be estimated by the filter, with 
which the non-Markovian dynamics of the qubit system can 
be obtained by tracing out the ancillary systems. 

The remaining contents are organized as follows. We 
briefly review the model of Markovian quantum systems in 
Section [H] In Section [Hll multi-ancillary systems driven by 
white noise are introduced for generating multi-Lorentzian 
noise. In Section lTVl a description of the qubit system is given. 
In Section [Vl we show that the qubit system satisfies a QSIDL 
disturbed by the colored noise from the noise model which is 
a combination of Lorentzian noises. A quantum Alter for the 
non-Markovian qubit system is discussed in Section |VT] An 



illustrative example is given in Section IVIII Conclusions are 
drawn in Section IVIIII 

II. Brief review of Markovian quantum systems 

A. White noise field 

A Markovian quantum system refers to a quantum system 
interacting with white noise fields. The white noise field can 
be dehned as 

h(t) = —= / &(w)e (1) 

satisfying the delta commutation relations 

[b{t),b\t')] =S{t-t'), [b{t),b{t')] = 0, (2) 

where the operator b is an annihilation operator of the field on 
the Fock space This white noise field may be described 
as a quantum stochastic process. With the dehnition O, 
an integrated operator can be dehned as Bt = b(t')dt' 
whose adjoint is bJ = b'^(t')dt'. They satisfy [Bt,Bl] = 

mm(t,t'), \Bt,Bt'] = 0 and thus the operator Qt = Bt + b\ 
is the quantum analog of the Wiener process and q{t) — 
h{t)+h^{t) is quantum white noise. Note that we have assumed 
that the initial state of the held on the Fock space S' is a 
vacuum state such that this process is analogous to Gaussian 
white noise with zero mean. 

B. Dynamical equation of the Markovian quantum system 

Considering a quantum system interacting with the white 
noise held, a QSDE for an arbitrary operator X of the 
quantum system can be written down to describe its Markovian 
dynamics as 

dXt = {-][Xt,Hs{t)]+CUXt))dt 

+dBl[XuLt] + [L\,Xt]dBt ( 3 ) 

with a generator 

g{X) = -\[X,Hs]+CL{X), (4) 

where L is the coupling operator of the system and £.(•) 
dehnes a Lindblad superoperator which can be calculated as 
CM {O) = iiV'i' \0, iV] + ^ [W't', 0]N for two arbitrary operators 
N and O with suitable dimensions. The two terms in the 
hrst row on r.h.s. of Eq. © describe the free evolution and 
the dissipation process, respectively. And the terms in the 
second row describe the inhuence of the white noise held on 
the system. Such an equation describes the dynamics of the 
system driven by an external white noise held, which has been 
widely used in the analysis and control of Markovian quantum 
systems GU . Note that throughout the paper we assume h = 1. 

C. Input-output relations 

The input-output relation of the Markovian quantum sys¬ 
tem is an important issue for observing the dynamics of the 
system, which can be written as 

dBout(f) = Ltdt dBt. (5) 

This relation shows the output held dBout{t) not only car¬ 
ries information of the system but also is affected by noise 

dBt CD. 


D. (5", L, H) description 

To concisely describe the interconnection of Markovian 
subsystems, the (S, L, H) description of quantum systems ifTSl 
has been developed. In terms of the {S, L, H) description, the 
Markovian system introduced in the above subsection can be 
systematically denoted as 

G = iS,L,H), (6) 

where a scattering matrix S describes the input-output relation 
of helds passing through beam splitters, the operator vector L 
is a collection of system operators interacting with the external 
helds, and H is the system Hamiltonian. 

The {S, L, H) description can also concisely describe the 
interconnection among subsystems. When we consider the case 
that the input held of the second subsystem G 2 is the output 
held of the hrst subsystem Gi, they can be denoted as a series 
product G 2 <1 Gi. In addition, the case that two subsystems 
are assembled together without any other connections can be 
described by a concatenation product GifflG 2 ifTSl . With these 
basic notations, a quantum feedback network can be described 
by an {S,L,H) description. 

E. Master equation 

The dynamics of the Markovian quantum system can also 
be described in the Schrodinger picture by using a master 
equation for the density matrix of the system pi in a Lindblad 
form as 

pl = -i[Hs,pl]+LUpl). (7) 

where the superoperator £*(•) is the adjoint of the Lindblad 
superopertor calculated as Clq{0) = ^A^[0, W^']-f O] 
for operators N and O with suitable dimensions. Note that we 
have assumed the held is in a vacuum state. And this Lindblad 
form master equation is a differential equation where the state 
variation only depends on the present state showing Markovian 
nature of the dynamics. 

III. Multi-ANCILLARY systems driven by white 
NOISE 

A. Dynamics of multi-ancillary systems driven by white noise 

To generate colored noise with a multi-Lorentzian spec¬ 
trum, we consider n ancillary systems driven by white noise 
in this section. We assume that the A:-th ancillary system is a 
Markovian linear quantum system which is described by an 
{S,L,H) description as 

G^ = (l,y/yllak,ujkalak), ( 8 ) 

e.g., an optical mode in a cavity, where ujk is the angular 
frequency and ak (a^) is the annihilation (creation) operator of 
the k-th ancillary system. Here the coupling operator is chosen 
as y/^Qk, where is a damping rate with respect to the 
white noise held ©I. Eor each ancillary system, the scattering 
matrix is an identity matrix, which means no scattering process 
for the helds is involved. Supposing the A:-th ancillary system 
is dehned on a Hilbert space f)§, it evolves on a Hilbert space 





The multi-ancillary systems can be defined as 

_ m m m m m 

a — bd O-Q bd • • • bb (_z^ bb ■ ■ ■ bb (_z^ 

= {1,M,Ha), fc = l,2,..-,n (9) 

on the space f)®" ® 5 ", where M = FA is a coupling 
operator for the multi-ancillary systems with respect to white 
noise with a damping matrix F = diag[Y^, • • • , 
and a collection of the annihilation operators for the multi- 
ancillary systems A = [oi, 02, • • • , a„]^. The internal Hamil¬ 
tonian of the multi-ancillary systems Ha can be expressed as 
Ha = A^VlA, where H = diag[a;i,a;2, • •' with angular 
frequency uJk for the fc-th ancillary system, k = 1,2 ,-•• ,n. 
Note that we assume that all the multi-ancillary systems are 
driven by the same white noise field. 

The unitary evolution for the multi-ancillary systems can 
be described by an evolution operator 0t in the interaction 
picture with respect to the white noise field satisfying a QSDE 
as follows 

d0t = {-{iHA+ ^M'iM)dt + dBlM-M'<dBt}Qt, (10) 

where dBt = [dB (t), dB (t), • • • , dB{t)\^ describes the white 
noise process. The generator for the multi-ancillary systems is 
GaiXa) = -i[Xa,HA] +CM{Xa), where Xa is an operator 
of the ancillary systems. 

Hence, a QSDE for the annihilation operators vector A for 
the multi-ancillary systems can be written as 

dA{t) = -{^ + in)A{t)dt-rdBt (11) 

with A{t) = 0 |A 0 t, where the variation of the operator A is 
driven by the white noise process dBt. 


We define 


C{t) = - 


FU(f) 


( 12 ) 


as a fictitious output. Then the fictitious output C{t) satisfies 
a QSDE as follows 

Ftp pip 

= -(—+i!^)^^(i)df+—dBt (13) 

whose formal solution can be expressed as 

Jto 2 

(14) 

with an initial state C{to). 


B. Multi-Lorentzian spectrum 

We have assumed that the multi-ancillary systems are a 
part of the environment such that the dynamics of the multi- 
ancillary systems are assumed to start from a long time ago so 
as to let to —>■ —00. Hence, a stationary version of C{t) can 
be obtained as 

C{t) = elCQt = f s(f - T)b(T)dT, (15) 

J —CO 

which is a convolution involving the white noise 
field b(t) = ,b{t)]'^ and a kernel 5 (f) = 


diag[Ci(f),6(f),--- with ^k{t) = 

fc = 1, 2, • • • ,n. 

The power spectral density for each component of the 
fictitious output C{t) is Lorentzian and calculated to be 

li 

S'fe(a;) =-2 - - - , k = ,n (16) 

-f + {to - UJkY 

where the center frequency Wfc and the linewidth 7^ are de¬ 
termined by the fc-th ancillary system’s the angular frequency 
and the damping rate with respect to the white noise field, 
respectively. The commutation relation for C{t) is determined 
by a memory kernel function, i.e., the Fourier transform of 
the spectrum, [C(f), C'l(f')] = J^~^[diag[ 5 'i(a;), • • • , 5 '„(u;)]], 
which is different from that of white noise. 


IV. Principal single qubit system 


A single qubit system is a basic unit of quantum compu¬ 
tation and quantum information, which is defined on a two- 
dimensional complex Hilbert space 1 ). The quantum informa- 


qubit which are denoted as |0) = 
respectively. 


1 excited states of a single 

0 

1 

and 1) = 

1 

0 J’ 


More generally, a density matrix p'^ is introduced to 
describe the state of an open single qubit system, i.e., one 
qubit system interacting with external environments or other 
quantum systems, which can be expanded as 

p'‘=+xax + yuy + z(Jz) (17) 

where 

0 1 ■ 

1 0 J ’ 

0 -i ■ 

i 0 J ’ 

1 0 
0 -1 


tXx - 


Cfz = 


are Pauli matrices and I is the 2 x 2 identity matrix and 
\x,y,zY' is the Bloch vector. For more details, see IT]. 

In addition, the ladder operators for the qubit system 

0 0 ■ 

1 0 J ’ 

0 1 ■ 

0 0 

are utilized to describe a state flip between |0) and |1), e.g., 
cr_|l) = |0) and (T+|0) = | 1 ). The ladder operators are also 
used to describe the interaction with external systems, e.g., in 
the Jaynes-Cummings model na. 

The Hamiltonian of the single qubit system we considered 
is given as 

Hs = (18) 

where ujg is the qubit working frequency. 


(T_ = 

cr+ = 
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Fig. I. Schematic diagram for the direct coupling between the multi-ancillary 
systems and the qubit system. 


V. Single qubit system interacting with 

MULTI-ANCILLARY SYSTEMS 

A. Dynamics of the augmented system 

In this section, we consider a general case that the single 
qubit system is strongly coupled with n ancillary systems via 
their direct interaction as shown in Fig. [T] where the multi- 
ancillary systems have dynamics as discussed in the Section 
m The augmented system is defined on an extended space 
f) 0 f)®" 0 T?. Note that the dynamics of the ancillary systems 
cannot be eliminated via the adiabatic elimination which is 
valid for the off-resonant case, i.e., there exists large detuning 
frequencies between the qubit system and the multi-ancillary 
systems ini. 

We assume that the interaction Hamiltonian for the cou¬ 
pling between the qubit system and the multi-ancillary systems 
is 

iF/= i(C'^E-S^C'), (19) 

where the direct coupling operator of the qubit system S can 
be expressed as E = [y/niai, y/K 2 a 2 ,- " , with 

the coupling strengthes y/Kk and the qubit system operators 
(Tfc, fc = 1, 2, • • • ,n. Note that C = —^ as given in (fTSl i. 

This augmented system can be described by using an 
{S,L,H) description as 

Gq-a = (I, M, Hs + Hi + Ha). (20) 

The evolution operator Ut of the total system satisfies a 
QSDE as follows 

dUt = { — i{Hs + Hi + HA)dt —-M^ Mdt + 

dBlM - M'<dBt}Ut. (21) 

Let X' denote any operator for the augmented qubit and an¬ 
cillary system. Its evolution can be defined as X). = U}x'Ut, 
which satisfies a QSDE written as 

dX' = -i[X'„Ht]dt + CM^i-^iXi)dt 
+ i[XfC}tt] + [tlCt,Xi])dt 
+ (dBl[X; Mt] + [Ml X,']dBt), (22) 

with Ht = Ul{Hs + HA)Ut, Mt = u}MUu Ct = ujCUu 
and Et = [//EC/*. 


In particular, for X' = X a qubit system operator, Eq. (l22l i 
reduces to 

dXt = -i[Xt, Hsit)]dt + Ei] + [El, Xt]Ct)dt (23) 

with Xt = UjxUt and Hs{t) = U^HsUf When X' = C, 
i.e., for the operator vector of the multi-ancillary systems, we 
have 

F'tT - riT - F'l'r 
dCt = -(— + m)Ctdt + —Stdf + —dBt. (24) 

Thus the solution of Ct can be written as 

Git) = Git) + l f 5(< - T)trdT. (25) 

^ Jto 

which shows the ancillary systems not only depends on the 
multi-Lorentzian noise vector C{t) but also is disturbed by 
the qubit system as indicated by the integral term in ( l25T l. This 
back action from the system to the ancillary systems will not 
happen in Markovian systems. 


B. Interaction picture with respect to the multi-ancillary sys¬ 
tems 

To show that the qubit system is driven by colored noise, 
we can move to the interaction picture with respect to the 
multi-ancillary system by defining an evolution operator as 
Vt = QlUt, whose evolution satisfies 

Vt = {-iHs- (EtC(f) - Clt)^)}Vt. (26) 

In this interaction picture, the system is described by 

Gl_, = (-, Hs + i(Ct(f)E - EtC(f))) (27) 

where C{t) has a Lorentzian spectrum. Hence, it is clearly 
seen that the system is driven by multi-Lorentzian noise in the 
interaction picture. 

Note that in the interaction picture, the system is driven 
by multi-Lorentzian noise as given in Eqs. (l26l l and dZTl l. 
The evolution of an operator X for the qubit system in 
the interaction picture is equivalent to that in the augmented 
system due to 

Vt^XVt = U}QtXQ\Ut = ulXQSlUt = u}xUt. (28) 

Hence, the operator evolution for the qubit system in Eq. ( |2^ 
is disturbed by multi-Lorentzian noise as well. 


C. Non-Markovian dynamics of the qubit system and its 
Markovian limit 

Substituting the solution (l25l l into (l23t . a non-Markovian 
Langevin equation for the qubit system can be obtained as 

Xi = -i[Xt,Hs{t)\+Clt)[XtM + pt,Xt]C{t) 

+79(5% Et)4Xt, Ei] + [El, Xt]79(5, E)* (29) 

where the convolution terms are expressed as 

I I 

D{S,^)t = - / S(f-r)E,dr. (30) 

Jto 
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Fig. 2. Schematic diagram for probing a non-Markovian quantum system. 

In particular, when (Ti = cr2 = • • • = (Tn = cr, i.e., every 
ancillary systems couples with the qubit system via an identical 
operator a, the non-Markovian equation (l29b is simplified as 

Xt = -i[Xt,Hsit)]+c^(t)[Xt,^t] + [d-l,Xt]c{t) 

+DiC,^^)t[Xt, at] + [alXt]DiC, a)t, (31) 


system as Ut which satisfies a QSDE as follows 

dC/i = {-i{Hs+Hi + HA)dt-]^M''Mdt 

-^L'fLdt + dB|M - M+dBt 

+dBlL- L'^dBt}Ut. (36) 

Then a QSDE for an operator X' of the augmented qubit and 
ancillary systems defined on () 0 1)®" can be derived as 

dX' = -i[X', Ht]dt + i[Xi C/St] + [ElCt,X't])dt 

+ + Ci^{X^))dt 

+ (dBt[x;,M(f)] + [Mt(Q,X']dBt) 
+d&t[KLt] + [LlX^]dBt, (37) 

where X' = Ujx'Ut, Ht = U^Hg + Ha)Uu Ct = U}CUt, 
tt = it = UlLUt, M{t) = ulMUt and dBt is the 

probing field process. 


where at = UjaUt- The kernel function (^{t) can be expressed 
as 

n 

(32) 

k=l 

whose corresponding power spectrum density 

n n 

='^'^kSk{u}) = ^-:^—- (33) 

is a combination of Lorentzian spectrum with weights Kk^k = 
1 , 2, • • • ,n determined by the coupling strengthes since ^k{t) 
has a Lorentzian spectra as in (fTfil l. The colored noise term 
c{t) can be expressed as 

^ V^^k{t - t)) b{T)dT, (34) 

which is driven by the white noise b{t). 

This Langevin equation OTT l coincides with the existing 
non-Markovian Langevin equations whose integral terms rep¬ 
resent the memory effect El, 11. Note that by choosing the 
parameters of the multi-ancillary systems, e.g., the angular 
frequency ojk and the damping rate jk, or the coupling strength 
Kk, the resulting noise spectrum ( |3^ can be approximately 
shaped as an arbitrary noise spectrum El- 

VI. Quantum filtering for non-Markovian 

QUANTUM SYSTEM 

A. The augmented system under a probing field 

To estimate the dynamics of the non-Markovian qubit 
system, a quantum filter can be constructed by using a probing 
field defined on a Lock space as shown in Eig. |2] The total 
system Gt can be described as 

Gt = ,Hs + Hi + Ha) (35) 

where L is the coupling operator of the qubit system for the 
probing field. We denote the evolution operator of the total 


Note that supposing an operator of the augmented system 
can be denoted as X' = Xq^Xa, the generator can be written 
as 

GriX') = gqiXq)(i^Xa+Xq(i^ga{Xa) 

(38) 

where 

GqiXq) = -i[Xq,Hs]+CL{Xq) (39) 

GaiXa) = -i[Xa,HA]+CMiXa) (40) 

are the generators for the qubit system and the ancillary 
system, respectively. 

One can write down the Langevin equations for the oper¬ 
ators of the qubit system. However, due to the commutation 
relations for the operators of the qubit system, these equations 
are nonlinear. Hence, it would be better to describe the 
augmented system by using a master equation. 

B. Unconditional Master equation 

By using the fact that the expectation of an operator X' 
in the Heisenberg picture is equal to that in the Schrodinger 
picture, we can obtain an unconditional master equation for 
the augmented qubit and multi-ancillary systems as 

Pt = —^[Hs + HA,pt]+C*M{pt)+£g{pt) 

+ [C^£,Pt] + [pt,£^G], (41) 

where pt is the unconditional state of the augmented system 
and the superoperator £*(•) is the adjoint of the Lindblad 
superoperator. 

As can be seen from Eq. (SD, the state evolution of the 
augmented system is Markovian, where the state variation 
only depends on the present state. One can also obtain the 
unconditional state pi of the qubit system by calculating 

Pt=Ha[pt], (42) 

which will not satisfy a Markovian evolution. Note that tr^ [•] 
means the partial trace with respect to the multi-ancillary 
systems. 




























C. Belavkin quantum filter 

Using the probing field, the system can be continuously 
monitored via homodyne detection, where a quadrature of the 
probing field is detected and can be used as an input to a 
quantum filter. 

It is easy to check that the probing field in a vacuum 
state satisfies a non-demolition condition which means the 
continuous measurement of the field does not change the 
observable of the qubit system m. Also, we assume the 
detection efficiency of the homodyne detector is perfect with 
100 % detection efficiency. 

Hence, we can follow an orthogonal projection approach to 
obtain a Belavkin quantum filter lfT9l . Il20l for the augmented 
system as 

diTtiX') = Trt{gT{X'))dt-iTrtiX'L + L^X')-7rtiX') 
XTTt{L + L'^)){dYt - TtfiL + L'^)dt) (43) 

where X' is an operator of the augmented system. Yt is the 
output field and dW = dYt — 'Kt{L + L'^)dt where W is called 
the innovation process and is equivalent to a classical Wiener 
process. The estimate of an observable X'^ is defined by a 
conditional expectation as X'^ = TTt{X') = E[X(| 3 ^t], where 
yt is a commutative subspace of operators generated by the 
measurement results Y{t), 0 < t < t. Note that the increment 
dW is independent of TTr(X'), 0 < t < t. 


D. Stochastic master equation 

The conditional expectation Trt(X') is defined for the 
augmented system and thus a conditional density pt for the 
augmented system can be defined by 

'Kt{X')=tY[ptX']. (44) 

Hence, a stochastic master equation for the augmented system 

can be obtained from the quantum filter as 

dpt = C/T(pt)df + (45) 

with 

Xl{M = Lpt+ptL'^-tT:[(LY Li^)pt]pt, (46) 

which is a Markovian stochastic master equation. The super¬ 
operator C/y is the adjoint of Qt- 

However, a stochastic master equation for the density 
operator pi of the qubit system is not in a Markovian form. 
Instead, we can trace out the ancillary system to obtain the 
conditional state of the qubit system pi as 


p1 = iT:a[pt]. (47) 

In practice, one cannot obtain an exact description of pi due 
to the infinite dimensional nature of the ancillary systems. 
However, a truncation can be made for the ancillary system, 
i.e, we can assume it is a A^-level system and thus it is possible 
to calculate an approximation to the partial trace (|47] |. 
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Fig. 4. The dynamics of observables for the qubit system in both non- 
Markovian and Markovian cases. The unconditional and averaged conditional 
expectation of observables in the non-Markovian case are denoted as (c. (i)) 
(green lines) and (blue lines), respectively. For the Markovian qubit, 

the unconditional expectation of observables are denoted as (cr"^{t)) (red 
lines). 


VII. An illustrative example 

In this section, an example of a single qubit system coupled 
with one ancillary system (i.e., n = 1) which converts white 
noise to Lorentzian noise is given in Fig. [3 Here, the direct 
and field coupling operators are specified as E = yTcTcry and 
L = y/fiqax, respectively. The corresponding parameters are 
set as uji = ujq = lOGHz, ki = 1, and 7 ^ = 0.8. The damping 
rate of the ancillary system with respect to the white noise field 
is 7 i = 0.6. The single qubit system is initialized in a state 

5(1 + O'x)- 

Fig. m shows the evolution of both the unconditional and 
averaged conditional expectation values of the observables ax, 
ay, and cr^ for the qubit system. The conditional state pt for 
the augmented system can be obtained from the quantum filter 
(l45l l and thus the conditional expectation of observables for 
the qubit system can be calculated as {a) = ® I)pt], 

where a is an observable of the single qubit system, e.g., a 
can be ax, ay or a^. Here, I is the identity matrix defined on 
the Hilbert space of the multi-ancillary systems. The averaged 
conditional expectations {dx,y,z) are plotted as blue lines. 

































which are obtained from the average for 500 realizations of 
the trajectories. The green lines represent the unconditinal 
expectations {ax,y,z) which are obtained from the results of 
the master equation of the augmented system dTlI i. It shows 
the quantum filter can estimate the non-Markovian evolution 
of the single qubit system. 

Compared with the non-Markovian trajectories, the uncon¬ 
ditional expectation values of the observables dx, <Jy, and Uz 
for the qubit system in the Markovian case are also plotted 
as the red lines in Fig |4] where the qubit is directly open 
to the white noise field and the probing field. In this case, 
the system dynamics obeys a Markovian master equation as 
pI = -'AHs,pV\ + It shows that not only 

the qubit in the Markovian case damps faster than that in the 
non-Markovian case but also the stationary states of the qubit 
in the two cases are different. 

VIII. Conclusion 

In this paper, we have investigated a non-Markovian quan¬ 
tum system in an extended Markovian representation frame¬ 
work, where multi-ancillary systems are introduced to convert 
white noise to colored noise with multi-Lorentzian spectrum. 
The multi-ancillary systems of this model play the role of the 
internal modes of the environment resulting in non-Markovian 
dynamics of the qubit system. Such a model is also compatible 
with methods of quantum control theory so that a quantum 
filter can be constructed to estimate the state of the non- 
Markovian system. An illustrative example involving qubit- 
cavity systems has shown the quantum filter can estimate the 
non-Markovian dynamics of the qubit system. In principle, 
our multi-ancillary model can approximately capture a non- 
Markovian environment with an arbitrary spectrum by redis¬ 
tributing the multi-Lorentzian spectrum. Then, a robust quan¬ 
tum filter for any non-Markovian systems can be constructed 
in future work. 
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